Dynamics of Josephson junctions and single-fiux-quantum networks with 
superconductor-insulator-normal metal junction shunts 
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Within the framework of the microscopic model of tunneling, we modelled the behavior of the 
Josephson junction shunted by the Superconductor-Insulator-Normal metal (SIN) tunnel junction. 
We found that the electromagnetic impedance of the SIN junction yields both the frequency- 
dependent damping and dynamic reactance which leads to an increase in the effective capacitance 
of the circuit. We calculated the dc I-V curves and transient characteristics of these circuits and 
explained their quantitative differences to the curves obtained within the resistively shunted junc- 
tion model. The correct operation of the basic single-flux-quanta circuits with such SIN-shunted 
junctions, i.e. the Josephson transmission line and the toggle flip-flop, have also been modelled. 

PACS numbers: 74.50. +r, 84.30.-r, 73.40. Gk 



INTRODUCTION 

Recently, the Josephson digital circuits operating on 
Single Flux Quantum (SFQ) pulses 0,0 have been con- 
sidered as electronic circuits suitable for integrating with 
Josephson qubits (see, for example, Refs. 0, 0)- Ap- 
plying Josephson SFQ electronics for the control and 
readout of the qubits has many advantages, including 
a high speed of operation and low operating tempera- 
ture in combination with rather small dissipating power, 
allowing a sufficiently close location of the elements on 
the chip. Moreover, this Rapid SFQ (RSFQ) electronics 
makes it possible to process the input and output signals 
of the qubits directly on chip. This possibility can extend 
the class of algorithms to be realized in the Josephson 
quantum computer. Finally, a similar Josephson junc- 
tion fabrication technology for both, the RSFQ circuits 
and the qubits, is an essential prerequisite to reach a joint 
circuit architecture. 

The operating principle of the RSFQ circuits is based 
on the single 27r-leaps of the Josephson phase tp in the 
overdamped junctions. Due to large damping, the driven 
Josephson junctions never switch completely into the 
phase-running regime with a large average voltage across 
the junction, but generate short SFQ pulses, V(t) = 
($o/2ir)d(p/dt, with quantized area, J V(t)dt — 5>o, 
where 3>o = h/2e « 2.07 x 10~ 15 Wb is the flux quan- 
tum [lj. For the most manufacturable Superconductor- 
Insulator-Superconductor (SIS) tunnel junctions with a 
not very large critical current density j c , operating at a 
temperature which is notably below the critical tempera- 
ture of the superconductor, (T < 0.5 T c ), intrinsic damp- 
ing is very small. Sufficient damping is therefore achieved 
due to an external low-ohmic resistance R s shunting the 
junction, so the McCumber-Stewart parameter |5j,|y], 



(3 C = u?Ju? v = (2ir/$ )I c R 2 s C, 



(1) 



the plasma resonance frequency of the Josephson junc- 
tion, ujp = (27r/ c /<i>oC) 1 / 2 , is just slightly lower than the 
Josephson characteristic frequency uj c — (27r/<I>o)/ c i?s 
and the plasma oscillations are strongly damped. 

Contrary to this, the operation of Josephson qubits re- 
quires vanishing damping in both the qubit's junctions 
and in the environment, including the coupled control 
and readout circuits. The effect of unsuppressed damp- 
ing is, however, dramatic and appears in fast decoherence 
of the qubit. Especially, the noise resulting from damp- 
ing at frequencies around the characteristic frequency of 
the qubit, Q, = (E - Ei)/% (typically, about 10 GHz), 
causes intensive relaxation, while the low-frequency com- 
ponents of the noise lead to a dephasing of the qubit (see, 
e.g. the review jjj). At short decoherence times, qubit 
manipulation is unfeasible. This situation can to a cer- 
tain degree be softened by weakening and/or switching 
on and off the coupling [g and by operating the qubit 
in the optimal points where the qubit is immune to the 
external noise in the linear order 0, ^} . In the case of 
the RSFQ circuit interface coupled to the qubit, the low 
resistance of the junction shunts R is the source of large 
broad-band current noise (oc i? _1 ) acting on the qubit. 
Moreover, this noise is generated by the resistors even in 
the quiescent (zero-voltage) state of the Josephson junc- 
tions. So, the problem of reducing the noise of RSFQ 
circuits coupled to the qubit has to be solved radically. 
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where I c is the critical current and C is the junction 
capacitance, is small, i.e. (3 C 2. At such values, 



FIG. 1: Electrical circuit diagram of (a) an resistively shunted 
and (b) SIN-shunted SIS Josephson junction (JJ). 
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A promising approach based on the possibility of 
frequency-dependent damping of SIS junctions by means 
of their shunting with Superconductor-Insulator-Normal 
metal (SIN) tunnel elements (see Fig. 1) has recently been 
proposed in Ref. (ll) . In particular, it was shown that 
shunting by a non-linear resistance drastically improves 
the noise characteristics of the circuit and ensure suffi- 
cient damping to achieve an almost non-hysteretic shape 
of the I-V curves. The applied model was, however, 
rather simplified and the special aspects of dynamical 
processes not investigated. In this paper, we address the 
problem of a microscopic model describing the SIS+SIN 
circuit, as well as the analysis of the dynamics of the ba- 
sic RSFQ networks within the scope of this model. The 
aim of this work was to demonstrate the full functional- 
ity and applicability of these circuits for the Josephson 
qubits. 



by Prober et al. |l2j, who approximated the nonlinear 
conductance term by three straight line segments. 

Equation (J2J is, however, too rough and does not de- 
scribe correctly the processes as the decay of plasma os- 
cillations caused by a short kick of the phase in the zero- 
current-biased circuit. In fact, after fast exponential de- 
cay and approaching the sufficiently low level (ss V g ), 
the amplitude Va of the oscillations V(t) — Va(*) cos ui p t 
starts to decrease much slower, because of small damping- 
being available for the small instant values of voltage V 
according to Eq. On the other hand, at sufficiently 
small (3 C < 1, these oscillations should decay very fast, in- 
dependent of their amplitude. To avoid this discrepancy, 
a more elaborate (microscopic) model of the circuit in 
Fig. lb has to be applied. 

MICROSCOPIC MODEL 



SIMPLE MODEL 

Due to the very large values of the zero-bias differential 
resistance and the strongly nonlinear dc I-V character- 
istic of the SIN junctions at low temperature (T <C T c ), 
their noise is small in the range of frequency up to 
Lu g = A/h, where A is the energy gap of the super- 
conducting electrode of the SIN junction. If the qubit's 
characteristic frequency is sufficiently small, Q < u> g , the 
influence of this noise on the qubit is expected to be 
weak. On the other hand, if uj g is lower than the charac- 
teristic Josephson frequency ui c = (2tt /$o)I c R, where R 
is the asymptotic resistance of the SIN junction at large 
bias (V ^> V g = A/e), a large effective damping of the 
Josephson junction is still possible at sufficiently small 
values of parameter (3 C given by Eq. Q with replacement 
R s — > R jll| . The capacitance C in this equation is 
equal to the sum of the capacitances of the SIS and SIN 
elements, Csis + Csin- 

For an analysis of the I-V characteristics of the SIN- 
shunted Josephson junctions the simplified equation of 
motion with nonlinear conductance term was numerically 
solved in Ref. viz., 

C^+l£^(V)+I cS m^ = I, (2) 
where the instant voltage 

Here, ^si c n(^0 ^ s the dc I-V characteristic of the SIN 
junction with the instant voltage V(t) as an argument. 
Equation J5J is a modified Resistively Shunted Junction 
(RSJ) model equation 0, 0], relating the terms taken 
at a given instant. Earlier, a similar model was applied 
for describing the quasiparticle damping in SIS junctions 



We will derive the current through the SIN junction 
starting from the tunneling Hamiltonian of Cohen, Fali- 
cov and Phillips [l^ . 

H = H + H T = H s + H N + eVN s + H T . (4) 

Thereby, term H$ (Hn) is the Hamiltonian of the super- 
conductor (normal) electrode; term eVN$ includes the 
finite voltage V across the junction. The tunneling is 
described by the Hamiltonian 

H T = Y,T k A d «+ T k, q Ckd\ (5) 

k,q 

and is considered to be a small perturbation. Here Cfe 
and c\ are, respectively, the destruction and creation 
operators for an electron in state k in the supercon- 
ductor electrode, while d q and d q are the corresponding 
operators for the normal electrode. The nonzero (gen- 
erally, non-constant) voltage V leads to the additional 
ac phase factors for, say, the superconductor-electrode 
operators, Ck — > Cfce 1 *^')/ 2 an d c\ — > cj t e~ ilp W/ 2 , with 
<p(t) = (2n/<P )f t V(ti)dt>. 

The number operators are equal to 

n s = c i° k > n n=J2 d i d i> ( 6 ) 

fc q 

so the tunneling current is expressed via the expectation 
values of the operators Ns and Nn, 

Ism = {i) = e{N s } = -e(N N ). (7) 

Applying the perturbation theory assuming adiabatic 
turning on from the past of the interaction Ht we ob- 
tain the standard first-order result (see, for example, 
Ref. H3), 

hm = -^J e + °V-*{[I{t),HTtf)])odt!, (8) 
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where (...)o denotes averaging over the ensemble Hq. In- 
troducing the retarded Green's function 

K(t - t') = i(2e/h 2 )8(t - t') exp[+0(i' - t)] 

x T k , q TZ,^([ci(t)dS),c k ,(t')d q ,(t>n 0l (9) 

k,q,k r ,<j' 

where 9(t) is the Heaviside step function, we finally arrive 
at the convolution integral expression for the current 



^SIN = 



/ K(t - t') sin 

J — OC 



d£ (10) 



with function K playing the role of a memory kernel. One 
can see that the tunneling current depends both on the 
instant value of phase ip(t) and on its values in the past, 
fit'), t' < t, so Ea. (|10|l describes the causal physical 
process. 

Actually, the relations similar to Eqs. © and (|10fl yield 
the quasiparticle component of the current in the micro- 
scopic theor y o f the Josephson tunneling developed by 
Werthamer |l5| and Larkin and Ovchinnikov |16| and 
presented in the time domain by Harris |l7j . For our 
case of the SIN junction, the kernel K(t) can be easily 
found from its Fourier transform giving the well-known 
dc I-V curve of the junction. In fact, in the special case 
of a constant voltage bias V = Vq the phase ip runs lin- 
early, i.e. ip = u> v t + const, where lo v = (2-7T /<&q)Vq, and 
Eq. l|lUfl yields the dc current as a Fourier integral, 



4m(V) = -4mi-V) = 



K(T)sm{w v T/2)dT. (11) 



The BCS-theory-based expression for Igj c N is given by the 
integral over the states in the energy representation, 



1 

eR 



\E\[f(E)-f(E + eV)} 
(E 2 - A 2 )V2 



dE, (12) 



where the actual range of integration is \E\ > A and 
where f(E) = (1 + exp (E /knT)^ 1 is the Fermi func- 
tion (see, e.g., Ref. 18] ). Applying the reverse Fourier 
transform to Eq. (|llfl we obtain 



K(r) = 



7sin( $ ow/7t) sin urdr, r > 0. (13) 



In the case of zero temperature, T 
curve has the hyperbolic shape, 



0, when the I-V 



rdc 



} (V) = [29(V) 1}9(\V\ - V g )(V 2 - V 2 f 2 R-\ (14) 
the integral in Eq. (|13H can be computed explicitly yield- 



ing 



where 



K(t) = -{2fi/eR)5'{T) + Ki(j), (15) 



Ki(t) = -(A/eR)6(T)Ji(u g T)/T. (16) 
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FIG. 2: Real (solid line) and imaginary (dashed line) parts 
of the product = uj\Y' (lo) +iY"(ui)] calculated for the 

case of zero temperature, T = 0. The solid line curve also 
shows the shape of the dc I-V curve of the SIN junction (see 
the first equality in Eq. 11178 ^. The increase in temperature 
T leads to a rounding of the sharp corners of these curves at 
uj = ±uj g . The inset shows the equivalent electrical circuit of 
the SIN junction, which comprises (from left to right) the dy- 
namic bias-dependent conductance, the dynamic capacitance 
(both frequency dependent) and the geometrical capacitance 
of the junction barrier. 



Here 8' is the time derivative of the Dirac delta-function 
and Ji is the Bessel function of the first order. The 
first term in Eq. (|15l) describes the linear component of 
damping due to the ohmic asymptotic at high voltage, 
| V\ ^S> V g , while the second term, given by Eq. llflj l. 
describes the dispersive damping due to a strong non- 
linearity of the I-V characteristic in the vicinity of 
V g . (Compare the shape of the kernel Eq. (|16f) with 
(7rA 2 /TieR)9(T)Ji{<jjgT)Yi(tjjgT), where Yi is the Bessel 
function of the second kind, obtained by Harris ^} f° r 
the SIS junction.) In the case of nonzero temperature, 
T ^ 0, both integrals in Eqs. lfT2l and fHfl) yielding the 
kernel K(t) should be calculated numerically. 

The sine Fourier transform of the kernel K(t) Eq. (jl 1|) 
and the corresponding cosine transform give the real and 
imaginary parts of the junction's complex admittance, 
respectively, Y = Y' + iY" , viz., 



uY'(w) = -4 c N (^/2e) 



K(t) sin tordr, (17) 



K(t)(coslut - I) dr. (18) 



These relations arise from Eq. (|10J) on the assumption of 
a small ac voltage V — v a cos cot yielding ip = a sin tot 
with a = 27rw a /ci;$o <C 1- The corresponding real and 
imaginary parts of the admittance Eqs. I|17|) and i|18|) de- 
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termined by the casual kernel K(t) obey the Kramers- 
Kronig relations. 

As expected, the nonlinear dc I-V curve igDj(V) yields 
according to Eq. Ijl7(l the frequency-dependent damping 
Y'(ui). The odd function Y"(u) describing the junction's 
reactance is positive for u> > 0, so the SIN tunnel junc- 
tion behaves as an equivalent frequency-dependent dy- 
namic capacitance C(u>) = Y"(u>)/lo with the frequency- 
dependent losses Y'(ui), which is added up with the geo- 
metrical capacitance of the sandwich Csin- For the case 
of zero temperature, T = 0, the plots of the components 
Eqs. (fT7|) and (|T%)) are presented in Fig. 2, where the inset 
shows the equivalent electrical circuit of the SIN junction. 
Note that at u> — > 0, the value of the dynamic capaci- 
tance is finite, i.e., C = (2u> g R)~ 1 . For rather transpar- 
ent barriers (for example, for the SIN junctions with an 
Al superconducting electrode and a specific resistance of 
barrier p = 30 fl ■ /xm 2 ) this capacitance is comparable 
with the geometrical capacitance of the barrier (about 
50fF/>m 2 ). 

EQUATION OF MOTION AND ITS SOLVING 

For computing the SIS junction current in the circuit 
Fig. lb we can naturally apply either the microscopic 
tunnel model fl5t fl6| in the form of 01 or the sim- 
pler adiabatic model of the Josephson junction giving 
-feis = Ic shup 19]. In our case of significant total damp- 
ing due to sufficiently low resistance R of the SIN junction 
the quasiparticle current of the SIS junction can be 
neglected. Moreover, since the characteristic frequency 
lo c is appreciably lower than the gap frequency of the SIS 
junction ujg IS = 2Asis/ft (Asis is the energy gap of the 
SIS junction electrodes, which in the most favorable case 
should be much larger than A [23), one can neglect the 
dispersion of the supercurrent. This dispersion is essen- 
tial at frequencies u> « wf* s and manifests itself as the 
logarithmic Riedel peak [2lj . Therefore, the model com- 
bining the microscopic description of the SIN junction 
and the adiabatic description of the SIS junction (with 
constant amplitude of supercurrent I c ) is adequate for 
the circuit Fig. lb. So, the equation of motion takes the 
form 

+ f K.jt-t')^-^ dt' = I. (19) 

J— oo 2 

Here we decomposed the kernel K (r) into two parts in 
accordance with Eq. I|15l) , presenting the asymptotic con- 
stant contribution of damping separately, by the second 
term on the left-hand side. 

Generally, the solving of an integro-differential equa- 
tion with a slowly decaying kernel (see Eq. (|16[0 is a dif- 
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FIG. 3: Fitting of kernel Ki(t) Eq. JTHJ (solid line) by the 
Dirichlet series Eq. 121)1 with three terms (dashed line). The 
inessential discrepancy of the curves at small values of time r 
is related to the imperfection of the fit of the Fourier trans- 
form Eq. 1171 at high frequencies, u ^> u g . Fitting of this 
Fourier transform, i.e. the dc I-V curve Eq. 1111 . in the most 
critical range of frequencies u) ~ uj g is shown in the inset. The 
dotted line curves in both plots show the shape of the kernel 
and the I-V curve given by two terms in the series Eq. IL'Dli 
(coefficient B2 = 0) that roughly corresponds to the case of a 
finite temperature of T ~ 0.3 T c . 



ficult task, because at each time step t — » t+At of the nu- 
merical integration one has to take a convolution integral 
over the time interval from —00 to (t + At), that dramati- 
cally slows down such calculations. However, the calcula- 
tions can be significantly accelerated if one approximates 
the kernel by a finite Dirichlet series , because the ex- 
ponential shape of the kernel makes it possible to avoid 
the time-consuming direct integration at each step. In- 
stead, only small corrections to the convolution integrals 
are computed at each time step. This procedure was re- 
alized by Odintsov et al. [l^ for the kernels derived in 
the microscopic model of the SIS junction |15l llq. For 
our case of an SIN junction having frequency dispersion 
of a relatively simple shape (see Fig. 2), the Dirichlet se- 
ries can contain only a few terms, but still describe the 
dynamics adequately, 

A N 

n— 1 

here coefficients B n and p n are the complex numbers with 
Re(p„) < 0. 

Figure 3 shows the result of the kernel fitting by the 
Dirichlet series Eq. I|20() with N = 3 for the set of frequen- 
cies: Impi = 0.95 ujg, Imp2 = <-^>g and Imj>3 = 1.05 u> g . 
One can see that the approximating function captures 
well both the behavior of kernel -Ki(t) and its Fourier 
transform (shown in the inset). Interestingly, the two- 
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FIG. 4: The autonomous I-V curves calculated within the 
frame of the microscopic SIS+SIN model and the RSJ model. 
The values of the McCumber-Stewart parameter j3 c for the 
SIN-shunted junction are 0.3, 0.5, 1, 2, 5 and 100 (solid 
lines), while for the RSJ model the values are 1, 2 and 3 
(dashed lines), in sequence from the top curve to the bottom 
curve. The dash-dot arrow indicates the switching from the 
superconducting to the resistive state due to the current bias 
regime. The bottom curve calculated for large j3 c = 100 S> 1 
practically coincides with the I-V curve of the stand-alone 
SIN junction. The fragment of this curve is also shown in the 
inset of Fig. 3 by dotted line. The value of the characteristic 
voltage is equal to V c = 1.5A/e in this plot as well as in the 
plots which follow. 



term approximation of the series Eq. I|20|) with the fre- 
quencies Impi = 0.95 Lo g and Imp3 = 1.05 uj g (shown by 
dotted lines) yields a reasonable approximation of the 
dc I-V characteristic for finite temperature. The ob- 
tained set of coefficients B n and the damping factors 
Rep n in Eq. (|20[1 made it possible to reduce the numer- 
ical solving of the integro-differential equation Eq. 1|19|) 
to almost that of an ordinary differential equation. For 
the single junction circuits the simulations were tech- 
nically performed by applying the fourth-order Runge- 
Kutta method. For the multi-junction circuits the sim- 
ulations were done with the help of the fitted program 
code PSCAN developed earlier for both the RSJ and the 
tunnel junction models USES!- 

JOSEPHSON JUNCTION NETWORKS WITH 
SIN SHUNTS 

Applying the described procedure of solving the equa- 
tion of motion, we first found the dc I-V characteris- 
tics of the circuit Fig. lb biased by a constant current. 
The resulting curves are shown in Fig. 4, where they 
are compared with the I-V curves given by the RSJ 
model (Fig. la). One can see that the shapes of these 
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Time x = cot 
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FIG. 5: (a) Electric diagram of the SIN-shunted Joseph- 
son junction inserted in the superconducting loop with the 
dimensionless inductance Pl = 2 and driven by a step- 
pulse of the magnetic flux, (b) Switching characteristics of 
the SIN-shunted (solid lines) and resistively shunted (dashed 
lines) Josephson junctions induced by the short pulse (dot- 
ted line) of flux of unit amplitude (<& e (t)/&o = 0.1 at uj c t < 
10 and 1.1 at ui c t > 10) for two sets of values of /3 C and /3l- 
The kernel K\ is approximated by two terms in Eq. 12011 . 



curves are qualitatively similar, although the values of 
[3 C w 0.3 — 0.5 ensuring sufficiently small hysteresis in 
the curves for the SIN-shunted junction are appreciably 
smaller than the corresponding values in the RSJ model 
{[3 C ~ 1 — 2). Moreover, the former curves exhibit char- 
acteristic plateaus at V < A/e. Interestingly, the size of 
such plateaus developed at small (3 C is somewhat smaller 
than that obtained in the simplified (phcnomenological) 
model of the SIN junction (cf. Fig. 1 of Ref. HH). 

The dynamical process of the jump of the Josephson 
phase by 2tt was modelled in the circuit shown in Fig. 5a. 
The inductance of the superconducting loop L closing the 
shunted SIS junction is comparable with the Josephson 
inductance, i.e. the dimensionless parameter 

L = 2ttI c L/$ (21) 

is equal to 1 — 2. A short step-pulse of magnetic flux of 
the magnitude of <I>o was applied to the loop. The tran- 
sient behavior caused by this pulse is shown in Fig. 5b. 
One can see that for the case of sufficiently large damp- 
ing, p c — (5l — 1, the curves for the microscopic SIN and 
phenomenological RSJ models practically coincide. For 
smaller damping, i.e. for the values (3 C — 4, (3l — 2 giving 
the bare resonance frequency of the circuit Fig. la equal 
to luq = ujdp- 1 ^ + fil 1 )] 1 ' 2 w 0.61w c « 0.92w ff < uj g , 
the difference between the corresponding curves is sub- 
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FIG. 6: (a) Electric circuit diagram of the 4-cell Josephson 
transmission line constructed from SIN-shunted Josephson 
junctions with identical critical currents 7 C . The McCumber- 
Stewart parameter of these junctions is equal to j3 c = 0.3. 
The bias currents lb are equal to 0.9 J c . The value of each 
inductance is L = 0.4 <E>o/(27r7 c ). (b) The SFQ voltage pulses 
on the JTL cells, corresponding to the propagation of the 2n- 
leap of the phase; the voltages across different cells are offset 
for clarity. 



stantial. The transient behavior of the SIN-shunted junc- 
tion shows oscillations of larger amplitude and smaller 
frequency, compared to the RSJ model. The reason for 
this behavior is the effect of the dynamic capacitance 
C(uj) of the SIN junction (see Fig. 2), giving uj c RC « 
1.4. This leads to a 20% decrease in the resonance 
frequency ujq. The corresponding effective value of the 
McCumber-Stewart parameter f3 c becomes therefore re- 
spectively larger. 

We also present the results of simulations of the ba- 
sic RSFQ circuits, i.e. the Josephson transmission line 
(JTL) and the toggle flip-flop (TFF) Jl|. Figure 6a shows 
the electric circuit diagram of the JTL consisting of the 
chain of SIN-shunted Josephson junctions connected in 
parallel by relatively small superconducting inductances 
L. An SFQ step pulse of flux is applied to the leftmost 
loop of the line and causes sequential triggering of 27r- 
leaps in the junctions Jl, J2, J3 and J4. As a result, 
the SFQ voltage pulse is transferred along the line with 
a small time delay, ~ 2it/uj c , on each cell (see Fig. 6b). 

Figure 7a shows the electric circuit diagram of the 
TFF. Due to the appreciable value of the storing induc- 
tance L3, this circuit has two stable states that differ 
by the direction of the dc current circulating in the loop 
of the interferometer J1-L3-J4. The SFQ pulses arriv- 
ing at the input port cause an alternating switching of 
the TFF, while the auxiliary junctions J2 and J3 pre- 
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FIG. 7: (a) Electric circuit diagram of TFF constructed 
from the SIN-shunted Josephson junctions Jl, J2, J3 and 
J4 having the critical currents 7 c i = 1.35 I c o, 7 C 2 = 1.17 7 c o, 
I C 3 = 1-55 Ico and 7 C 4 = 1.4 7 c o, respectively. The values of 
inductances L0, LI, L2, L3, LA, L5 and Lg expressed in units 
L u = $o/(27t7 c0 ) are 0.4, 0.2, 0.2, 3.0, 0.8, 0.8, and 0.05, re- 
spectively. The bias current 76 is equal to 1.31 I c o- (b) The 
time dependence of the current flowing through the storing 
inductance L3 (solid line) and the 2-7r-leaps of the phases on 
the junctions Jl and J4 (dashed lines). Here the frequency 
ui c — (2ir /&o)I c oR, where the tunnel resistance 7? is similar 
for all SIN junctions. 



vent a back-reaction of the interferometer on the SFQ 
pulse source 0. The switching of the direction of the 
current circulating in the interferometer loop leads to a 
rectangular pulse of flux induced in inductance L2> (see 
Fig. 7b). Such pulse, as a control signal, can, for exam- 
ple, be applied to the loops of the Josephson flux qubit 
of the double SQUID configuration [25| . 

DISCUSSION 

In summary, we have developed the microscopic model 
of the Josephson junction shunted by the nonlinear el- 
ement based on the SIN tunnel junction. The behav- 
ior of the circuit is adequately described by the integro- 
differential equation with a memory kernel reflecting the 
casual dependence of the tunneling current on the phase 
across the junction. In contrast to the simplified model 
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proposed in Ref. , this model captures the expected 
features in the behavior of the SIN-shunted Josephson 
circuits, including a clear dependence of the decay of 
plasma oscillations on their frequency. 

Besides the frequency-dependent damping, the SIN- 
junction shunts have the effect of an enhanced capaci- 
tance. The resulting capacitance includes both the ge- 
ometrical capacitance of the junction and the dynamic 
frequency-dependent capacitance that affects the shape 
of the I-V curve increasing the hysteresis. That is why 
the calculated I-V curves are almost similar to those 
given by the RSJ model for somewhat larger values of 
the McCumber-Stewart parameter, i.e. /?^ SJ = 1-3. Still, 
the values C rs 0.3 ensuring a sufficiently small hystere- 
sis in SIN-shunted Josephson junctions and the function- 
ality of the multi-junction circuits seem to be feasible. 
The most challenging condition to be met is the real- 
ization of high-quality SIN junctions with a high trans- 
parency of the tunnel barrier, having a specific resistance 
p < 30 Q ■ /im 2 ■ The recent experiments with Nb- 
A10 X -A1 (at T > 1.4K) and Al-A10 x -Cu (at T < IK) 
junctions have shown that the ratio of the zero-bias resis- 
tance to the asymptotic tunnel resistance can achieve suf- 
ficiently large values, i.e. > 30-50 HI, |22l. (This behav- 
ior of the SIN junctions can be roughly described by the 
two-term approximation of the memory kernel Eq. J20J-) 
Using such SIN junctions as shunts for Nb SIS junctions 
can make it possible to significantly reduce the noise of 
the circuits in the quiescent state in the frequency range 
up to A.A\/h ~ 50 GHz, i.e. within the working frequency 
range of the Josephson qubits of different types. 

Finally, the modelled behavior of the simple RSFQ net- 
works (JTL and TFF) is qualitatively similar to that of 
the conventional (resistively shunted) circuits. A prelim- 
inary evaluation of the ranges of functionality was found 
to be quite good. The simulation of more complex cir- 
cuits can be easily performed using the slightly modified 
code PSCAN developed earlier for the tunnel junction 
model [HEl. So, RSFQ networks with SIN-shunts im- 
plemented in the shell that surrounds a Josephson qubit 
core offer a promising approach to achieve joint RSFQ- 
qubit operation. 
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